A FORMULA FOR A QUARTIC INTEGRAL: A SURVEY OF 
OLD PROOFS AND SOME NEW ONES 
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Abstract. We discuss several existing proofs of the value of a quartic integral 
and present a new proof that evolved from rational Landen transformations. 
We include our personal renditions as related to the history of these particular 
computations. 



1. Introduction 

The evaluation of definite integrals has attracted the scientific community, both 
professional and amateurs, for a long time. Some of these evaluations are collected 
in tables. The authors' favorite one is by I. Gradshteyn and I. Ryzhik [TT]. It was 
the naive attempt to prove all the formulas in [TT] that produced [6|. A recent 
review of this book [lOj describes how R. Feynman 'once claimed that he acquired 
his initial reputation not as a physicist, but as a redoubtable evaluator of integrals'. 

In the Fall of 1995, when the author (VHM) was teaching a graduate class in 
Analysis, a student (George Boros) stated that he could prove the identity 



(1.1) 



dx TT P-mifl) 



(* 4 + 2 ^ 2 + 1 )" l+1 2 [2(a+l)r+5 



where P m (a) is a polynomial in a. The integral above will be denoted by No^(a; m). 
This was the introduction of VHM to the wonderful world of definite integrals. 
George's original expression for the coefficients di(m) in 

m 

(1.2) P m (a) = Y,di{m)a l 

1=0 

was 

*«-SSl^(?)C , + *)("-"-; i )( , ;0C7-; i 

Clearly not very useful. Intrigued by the fact that symbolic languages, such as 
Mathematica or Maple, were unable to evaluate this integral, and being unaware of 
any useful technique available for the computation of integrals, the second author 
joined George Boros in the study of this polynomial. 
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Computing some of these coefficients we observed that they are all positive. For 
example, if m — 5, we obtain 



(1.3) {di(5) : < I < 5} 



J4389 8589 7161 777 693 63 
\ 256 ' 128 ' 64 ' "8"' ~16~' ~8~ 



In this paper, we describe several proofs of the identity (jl.ip where the coefficients 
di{m) are given by 

a., *M-*--t*ft:r)( m r)(f> 



2. A SURVEY OF PREVIOUS PROOFS 

2.1. The elementary proof. George's original idea was simple, but had profound 
consequences. The change of variables x — tan (9 yields 



7r/2 / 4/) \ m+l 



Nn Aa-.m) — l — -. ^ x 

' \ sin 4 9 + 2a sin 2 9 cos 2 9 + cos 4 9 



(W 



Now, we observe that the denominator of the trigonometric function in the inte- 
grand, is a polynomial in u — 29. In detail, 

sin 4 9 + 2a sin 2 9 cos 2 9 + cos 4 9 = (1 + a) + (1 - a) cos 2 u. 

In terms of the double-angle u = 29, the original integral becomes 

Ar I \ o-fm+l) I" ( (1 + C0SU) 2 \ du 

N Aa-.m) =2 [ - m+L > / 1 — x . 

' V ; 7 \{l + a) + (l~a)cos 2 uJ 1 + cosm 

Expanding the binomial (1 + cosu) 2m+1 , the reader will check that 

(2.1) / [(l + a) + (l-a)cos 2 up m+1) cos J udu = 0, 

Jo 

for j odd. The vanishing of half of the terms in the binomial expansion turned 
out to be a crucial property. The remaining integrals, those with j even, can be 
simplified by using the double-angle trick one more time. The result is 

N 04 (a- 1 m)=Y2-3[ I / [(3 + a) + (1 - a) cosi;]" (m+1) (1 + cosvY dv, 

\ 2j J J 

where v = 2u and we have used the symmetry of cosine about v — tt to reduce the 
integrals form [0,27r] to [0,7r]. The familiar change of variables z = tan(i;/2) pro- 
duces the complicated formula for the coefficients d/(m) given in the Introduction. 

Excited with this evaluation and considering that we could not find it in [IT] , 
we wrote this as a paper and submitted it for publication. The reviewer was un- 
forgiving: Anybody who knows anything about hyper geometric functions would be 
able to do this. We realized that it was time to learn about hypergeometric func- 
tions. Eventually, we convinced an undergraduate at Tulane University to check 
the details of the evaluation discussed above. They formed part of her final Honor 
Project and the results appeared in [7]. 
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2.2. The hypergeometric proof. The integral iVo,4(a; m) can not be found in 
the compendium [llj in its original form. Naturally, it is possible that it appears in 
some disguised form, after a change of variables. In this case, the simplest change 
is t = x 2 , that produces 

(2.2) N 0A (a;m) = - 



2 J {t 2 + 2at + 1)™+! ' 
Integrals of this type are listed in section 3.252 of [11]. Indeed, 3.252.11 states that 

poo f- 1 rff 

The special functions involved here are the gamma function defined by 

poo 

(2.3) T(s) = / f-Vdt, 



that the reader must have seen as a generalization of integer factorials (namely that 
r(n) = (n — 1)!), and the beta function defined by 

(2.4) B(x,y)= [ f-^l-ty^dt, 

Jo 

that appears in beginning courses in Statistics, and finally the associated Legendre 
function Pjf{z) which was not very familiar to us. The table [llj is a good first 
source for special functions; in particular, sections 8.7 and 8.8 are dedicated to 
Ptf(z), where one learns that these functions are solutions of the differential equa- 
tion 



y = 0. 



At the beginning of Section 8 in [IT], one finds the formula 8.702: 

Z + 1\ I 1 — Z 



(2-6) TO = rO^)^) 2 F l{ -v,» + l;l-»; 

The hypergeometric function 2-F1 (a, b; c; z) encountered above, is defined as the 
series 

(2.7) ^(a, bl c;z):=±^^z\ 
where 

(2.8) (a) A = a(a+l)(a + 2)---(o + fc-l) 

is the Pochhammer symbol. This function includes most of the elementary func- 
tions. For example, 

(2.9) (1 + z) n = 2 F 1 (-n 7 b;b;-z), for any b, 
and 

(2.10) \nz = {z- 1)2^1(1, 1; 2, l-z) 

just to name a couple. Section 9.121 of [llj contain many more. 
Using these identities in (jl.l[) produces 

N aA {a-m) = T n - 1 l 2 {a+l)-^ n+l ^B(2m+l,\) 2 F X [-m,m + l;m+§; 1=2] , 
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and the formula (|1.4|) follows from here. The details appear in [2]. The reviewer of 
[7] was correct: not much of the hypergeometric world is required to prove this. 

2.3. A detour into the world of Ramanujan. Not knowing about hypergeo- 
metric functions and before we search for this type of proof, we observed that 

f°° dx tt 1 

(2.11) 



o x 4 + 2ax 2 + 1 



o bx 4 + 2ax 2 + l 2V2^W5' 
From here one sees that if we define 

poo J 

(2.12) g(c) := ' 
and 

(2.13) h(c) = \Ja + Vl + c, 
then (|2.11[) shows that g(c) = TT\f2h![c). In particular, 

(2.14) h'(0) = J_JVb,4(o;0). 

7TV 2 

Now it is clear how to proceed: further differentiation yields 

(2.15) Ja + VTT^=Va~+l + ^J2^ N 4 (a; k - l)c k . 

n v 2 k=1 * 

Using Ramanujan Master's Theorem [I], we were able to establish the identity 



/ 4m\ / 2m 



m . 

\2m 



-\ -i 



x 7Vo,4(a; m). 



r°° ~.ro-l i 
(2.16) / . = i-2 6m + 3 / 2 

The integral on the left can be evaluated in elementary form to produce (|1.4p . The 
details appear in [4]. It was never clear to the second author, how George thought 
of using Ramanujan's Master Theorem. He simply knew how to integrate. 

2.4. Another elementary proof. The elementary proof described below was 
shown to us by an (anonymous) referee. Write (|2.14p as 

with 

(2.18) /(c) = ylb+WTTc-l). 
Differentiating m times with respect to c yields 

(2.19) m!(-l) m / — —— =7r%/2/ (m+1) (c). 

v ' v ; J (x 4 + 2(6 - l)x 2 + 1 + c) m+1 

Hence 

(2.20) iVo.4(fe-l;m) = ^(-l)"7 (m+1) (0). 

m! 

But 

1/2 



ly 
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Now use the expansion 

(-l) l (n + 2i- 1)! /cy" 



(^-.)-(i)-(. + -g t - i, r + r„-" , a) 

given in [TT], Formula 1.114, and the identity 



y ' \n) n2 2 "-! \n-l 

and taking the coefficient of c m+1 to obtain 

r 2m — 2i\ ( m + i 

I 

»=o 

as desired. 



(2.22) N 4 (a; m) = 7 T77 V ' f + ! 2>+l)', 



3. A CHANGE OF VARIABLES AND THE NEW PROOF 

It soon became clear that the change of variables 

x 2 -l 



(3.1) 



y = R 2 (x) 



2.T 



was the key to extend the elementary method to integrals of higher order. The 
inverse has two branches 



(3.2) 



where the plus sign is valid for x G [0, +oo) and the other one on (— oo,0]. The 
rational function R 2 arises from the identity 



(3.3) 



cot 29 = R 2 {cot0). 



This change of variables gives the proof of the next theorem and it leads to inter- 
esting equality among integrals. 

Theorem 3.1. Let / be a rational function and assume that the integral of / over 
R is finite. Then 



(3.4) 



f{x) dx 



f(y + vV + i) + f(y - vV + 1) 



dy 



f(y + vV + i) - f(y - vV + i) 



ydy 



Moreover, if / is an even rational function, the identity 
replacing the interval of integration from R to M + . 



remains valid by 



The reader will check that the integrand on the right hand side of (|3.4p is also a 
rational function. This new function is an even rational function if / is such. 

The first application of Theorem 13.11 was given in [3] where we proved that the 
integral 



(3.5) 



U 6 (a,b;c, d, e) 



cx 



f dx 2 + e 
c 4 + bx 2 + 1 



dx 
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is invariant under the change of parameters 
(3.6) a - 

b - 



ab-\- 


- 5a 


+ 56 + ! 


(a 


+ b 


+ 2) 4 /3 


a 


+ b 


+ 6 


(a-f 


b + 


-2) 2 /3 


c 


+ d 


+ e 


(a 4 


b + 


-2)2/3 


c(b- 


f2) 


+ 2d + 




a + b + \ 




c + 


e 


(a 4 


b + 


2)!/3' 



Iteration of (J3TBJ) produces a sequence of parameters (a„, b n ; c„, d n , e n ) with 

(3.7) U 6 (a n , b n ; c„, d n , e„) = U 6 (a, b; c, d, e). 

Moreover, as n — > oo, we have that a n , b n — > 3 and there exists a number L such 
that (c„, d n , e n ) — >• L(l, 2, 1). Passing to the limit in (13. 7|) we obtain 

7T 

(3.8) U 6 (a,b;c,d,e) = -L. 

This is precisely the rational function analog of the transformation of parameters 

(3.9) a + b 

that leaves the elliptic integral 

(3.10) G(a,b) = 



2 

ab. 



it/2 



dip 



\/ a 2 cos 2 f + b 2 sin 2 y> 

invariant. In this situation, the sequence (a n ,b n ) obtained by iterating (|3 . 9[) con- 
verges to a common limit: the arithmetic-geometric mean of a and 6, denoted by 
AGM(a, 6). In the limit, we have 

^ G ^ 6 ^2AGM(M)' 

precisely as in (13. 8|) . Information about the AGM can be found in [S] and [S]- 

This was the beginning of a series of results on changes of parameters that leave 
integrals of rational functions invariant. We have called (|3.6p a rational Landen 
transformation. These were extended to all even rational functions in [5] and the 
convergence of the iterations of these transformation were described in [13] . The 
extension to arbitrary rational functions on the whole line is described in [14] and 
|15j . It is an open question to produce rational Landen transformations for a general 
rational function on the interval [0, oo). The vanishing of the integrals such as (|2.1[) 
does not happen if the integrand is not an even rational function. 

From the point of view of numerical calculations, the rational Landen trans- 
formations give a procedure to evaluate definite integrals in iterative form. For 
example, in the simplest case of the elementary integral 

dx 

, ax 2 + bx + c ' 



(3.12) / = 
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with b 2 — Aac < 0, we have shown in [15j that I is invariant under the change of 
parameters 

(o„ + 3c„) 2 - ib 2 n 



(3.13) 



b n +i 

Cn+1 



(3a„ + c n )(a n + 3c„) - 



3(a 



6, 



(3a„ + c„)(a n + 3c„) - 5,2 

(3a» + c») 2 ~3b 2 
(3a n + c n )(a n + 3c 

1 ) 



Moreover, the error defined by 

(3.14) e„ = (a„ - \ \/ Aac - b 2 , b n , c„ - \\/ Aac - b 2 ), 

decays to at a cubic rate, that is, ||e n+ i|| < C||e n || 3 . 



The transformation ()3.1|) and its higher order analogues, given by the rational 
function R rn defined by 

(3.15) cot(m0) = R m {cot6), 

have produced interesting results for definite integrals of higher degree. See [H] for 
details. 

The second goal of this note is to use the Landen transformation (|3.ip to produce 
a new proof of the identity (jl.ljl and thus obtain the more explicit natural form of 
the coefficients di(m) given in (|1.4|) . 

Theorem 3.2. For m € N, let 



(3.16) 
Define 



Q(x) = 



(x A + lax 2 + 1)™+! ' 



Q(y + vV + i) + Q(y - vV + i) 
y 



Then 
(3.17) 
where 
(3.18) 



vV + i 
Qi(y) -- 



T m(y) = 



Q( y + vV+T) - Q(y - vV + i) 

T m (2y) 
2 m (l + a + 2y 2 ) m + 1 ' 

m + /u A 



fc=0 



m — k 



y 2k - 



Proof. Introduce the variable <fi = y + \/y 2 + 1. Then ?/ — yy + l = 
y — 4> ~ (j)^ 1 . Moreover, 

Qi(y) = [Q(4>) + QW 1 )] +|J T j(Q(0)-Q(r 1 )) 
= -^[^QW + W 1 )] 



and 
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The result of the theorem is equivalent to 

(3.19) 2 m (1 + a +!(</> - 0- 1 ) 2 )) m+1 S m (4>) = T m (0 - 

A direct simplification of the left hand side of (|3. 19|) shows that this identity is 
equivalent to proving 

(3.20) ^ = T m {<p - r 1 )- 

First proof. One simply checks that both sides of (|3.20p satisfy the second order 
recurrence 

(3.21) C m+2 - {<P 2 + 4>~ 2 )c m+ l + C m = 0, 

and the values for m = and m = 1 match. This is straight-forward for the 
expression on the left hand side, while the WZ-method settles the right hand side. 

□ 

Note. The WZ-method is an algorithm in Computational Algebra that, among 
other things, will produce for a hypergeometric/holonomic sum, such as (|3.18|) . a 
recurrence like (|3.21j) . The reader will find in p~7] and [l6j information about this 
algorithm. 

Second proof. In [12] . one finds the generating function 

k 

(3-22) B t {z) = $>*)fc-iTT> 

k>0 

where (a)& = a(a + 1) . . . (a + k — 1) is the Pochhammcr symbol. The special values 

(3.23) B_i(z) = , and B 2 (z) = ^ , 

are combined to produce the identity 



(3.24) - ? =L= {B^(z) n+1 (-z)" +1 S 2 (-z)" +1 ) = J2 (" k k 

Replace n by 2m and z by (Ay 2 )" 1 to produce 



2 V / T+F(2y) 2m 
We also have 



I U 2 ™+! +0-( 2m + 1 )) = J2 k] J 



k=0 



Z . 



Thus, 



fc=0 

v 2m t( 2m : k V 2k 

fe=0 



T m (4> — (j)^ 1 ) = T m (2y) 



k=0 



z 
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We conclude that 

(3.25) T m (^ - r 1 ) = , 1 = (^ 2m+1 + r (2m+1) ), 

and the result follows from (j> + cj)^ 1 = 2\Jy 2 + 1. 

We now prove The identity in Theorem 3.1 shows that 

(3.26) / Q(a)da:= / Qi(y)dy, 

Jo Jo 

and this last integral can be evaluated in elementary terms. Indeed, 

T m {2y)dy 



1 " / m + fc\ y°° (2y) 2fc dy 
2 m ^ o \m-k)J (l + a + 2y 2 ) m + 1 ' 

The change of variables y = ^^t" i gives 



co ,2k 



t 2k dt 



k=0 

and the elementary identity 



(l + t 2 ) m + i: 



00 t 2k dt n /2k\ (2m- 2k\ fm> ~ ! 



(1 + £2)m+l 2 2m + 1 \ kj\ m-k 



gives 



I Ql{y)dy ~2 2 ^[2(l + a)] m + 1 / 2 ^ o {m-k) 2 
This can be simplified further using 



k\ ok f2k\ /2m- 2k' 

k ) V m — k J \k 



-l 



(1- 



and the equality (|3.26|) to produce 

(3.28) f Q(y) = * ^ \r^Y ^ + ^ " f 1 (1 + «)*■ 

Jo 2 + [2(l + a)] m+1 /2 \ vTO yl vm -fc y ' v ; 

This completes the proof of (II. with the coefficients di(m) given in (|1.4| . 
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